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Abstract 

The tomographic probability distribution is used to decribe the kinetic equations for open 
quantum systems. Damped oscillator is studied. Purity parameter evolution for different 
damping regime is considered. 

1 Introduction 

There exist several formulations of quantum mechanics (see recent review in pQ). In all these 
formulations the state of a system is associated with mathematical structure (complex wave func- 
tion, density operator, path integral) which is different from standard probability distribution 
used in classical statistics. On the other hand there exists a possibility of specific representation 
of quantum mechanics used in recent works 013] where the quantum state is described by the 
standard probability distribution (tomogram or tomographic symbol of density operator). This 
representation of quantum states is called "probability representation". This representation can 
be used both for conservative and open quantum systems. In this work we focuse on the open 
quantum systems. We review the approach and clarify some technical methods to deal with the 
tomograms of observables (operators). The tomograms are generalized functions. In view of this 
the calculations with tomograms need some specific attention. 

We will study quantum kinetic equations which describe damping in quantum domain ( jUEl 
IH1IZ]). Evolution of conservative quantum system is described by Schrodinger equation for wave 
function |SJ . The quantum systems interacting with a heat bath (open systems) obey to quantum 
kinetic equations written for density operator introduced in OHO]. The density operator is the 
hermitian operator with nonnegative eigenvalues. The diagonal matrix elements of the density 
operator in arbitrary basis have the meaning of standard probability distribution. In view of 
this the trace of density operator has to be equal unity. Quantum kinetic equations must to 
be close in such form that the described general properties of the density operator have to be 
preserved in the process of quantum evolution. It means that the initial density operator p(0) 
is transformed into the density operator p(t) using a specific map. This map p(0) — > p(t) has 
to preserve nonnegativity, hermiticity and trace of the density operator. Such kind of map was 
introduced by Sudarshan ^T] for finite-dimensional Hilbert space and it was studied in [Ti]IT^ 
IT3] . The most general quantum kinetic equations which satisfy the demands to preserve the 
properties of density operator in the process of time evolution have been suggested in [bHITBI 



ITT)] . The evolution equation for damped harmonic oscillator was used in |17j . Recently [T^ITH] 
new probability representation of quantum mechanics was introduced. In this representation the 
quantum state is described by standard positive probability distribution (tomographic probability 
distribution). In the probability representation the Schrodinger (and von Newmann) evolution 
equation takes the classical-like form of equation for the standard probability distribution of Fokker- 
Planck type. On the other hand till now the general quantum kinetic equations for open systems 
have not been studied in the probability representation. The aim of our work is to obtain the 
kinetic equations preserving the properties of density operator in the process of time evolution in 
the probability representation of quantum mechanics. In fact we will use the specific procedure of 
tomographic star-product quantization j2fj] and geometric approach to tomographic map [21] which 
generalize the Moyal j22j star-product quantization. In the Moyal approach the quantum kinetic 
equation is written for Wigner quasidistribution function In the tomographic approach the 

kinetic equation is written for the standard positive probability distribution function. It is worthy 
to note that some nonlinear dynamical equations were considered recently [21] in the tomographic 
representation. The equations of quantum mechanics are linear equations though there are some 
attempts to introduce nonlinear equations describing quantum state evolution (see e.g. The 
linear generic quantum kinetic equations for open systems can be also presented in the form of 
Fokker Planck type equation as well. 

The paper is organized as follows. In next Section 2 we discuss some properties of tomograms as 
generalized functions. In Section 3 we review tomographic representation of quantum mechanics. 
In Section 4 we derive the tomographic form of quantum evolution equation. In Section 5 example 
of damped oscillator is given. Conclusion are presented in section 6. 

2 Tomographic symbol of the unity operator 

Let us consider the tomogram of unity operator. We check the validity of the obvious equality 

Wi * IVf = IOj. (1) 

This equality follows from operator equality 

i ■ i = i. (2) 

It can be shown that the tomographic symbol of the unity operator reads 

w t (X, fi, v) = -ir\X\5(ji)5(u). (3) 

The tomogram w^(X, fi,u) is a probability distribution if A is a density operator. For arbitrary 
operator A its tomogram can be negative and even complex function. We have to check the equality 
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First of all, we take the factors containing the variables X% and X 2 : 
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Taking into account the equality 
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we can reduce the equality under study to the following one 
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= -7i\X\5(/x)5(u). 

To evaluate this integral we make the change of variables defined by the formulas 

Mi + M2 = fyt, z/i + z/ 2 = fcz^, 



Mi - M2 = X, 



^1 — ^2 = A. 



The inverse relations read 



Mi = t^M + x), 
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M2 = ^(Ip-x), V2=^(kv-\). 



The Jacobian related to the introduced linear transformation has the form 
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It means that we get the desired result. 

We also check the validity of the equality 

for the tomogram of a pure state. This equality follows from the property 
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of density operator p^ of normalized pure state. The tomographic symbol w^(X, p,u) of the 
density operator p^ is described by the formula (see e.g. |20 ) 
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According to the formulas for star-product kernel given in [20] the star-product w$ * takes the 
form: 
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Choosing the factors containing the variables X% and X 2 and using the equality 
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Integrating the variables v\ and z/ 2 we get the expression 
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To obtain the result we make the change of variables defined by the formulas 

Hi + H2 = k, fn = ^{k + l), (20) 

fJLi-fJL2 = l, H2 = ^(k-l). (21) 



The Jacobian related to the change of variables reads 

d(/ii,/i 2 ) = 1 
d(k,l) 2' 

One can transform the expression for the star-product * to the following form 
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Taking into account the normalization of the wave function ip(x) 



we get the desired equality 
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If we take the arbitrary density operator p instead of the density operator p^ of pure state, 
we see that all formulas will be correct except the last one. Making corresponding substitution 

^(yiW(z 1 )i>(y 2 )^*(z 2 ) -> p(y 1 ,z 1 )p(y 2 ,z 2 ) we get 
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One can see that the following relation take place 

Thus, we proved the identity for the tomograms of pure quantum states. 
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3 Tomographic representation of quantum mechanics 

In this section we review the tomographic approach given in (THIGH!- In |2Z| an operator X = 
(Xi, X^r) is discussed for the case of N = 1 as a generic linear combination of the position and 
momentum operators X n = pnXn + Vnpn, where p n and u n are real parameters for n = 1, N, and 
X is Hermitian, hence observable. The physical meaning of /x = (pi, ...,Pn) and v = (ui, u N ) 
is that they describe an ensemble of rotated and scaled reference frames, in classical phase space, 
in which the position X may be measured. It was shown 27J that the quantum state of a system 
is completely determined if the classical probability distribution w(X, fx,u), for the variable X 
is given in an ensemble of reference frames in the classical phase space. Such a function, also 
known as the marginal distribution function or quantum tomogram, belongs to a broad class of 
distributions which are determined as the Fourier transform of a characteristic function (2H1- The 
following formula for a quantum tomogram was derived in |21j : 



w(X,n,v) = (5(X - nq - up)), 



(32) 



where (A) = Tr(pA) is the average value of the operator A on the state described by the density 
operator p. In j2H] it was shown that, whenever X is observable, w(X, /x, v) is indeed a probability 
distribution, as it is positive definite and satisfies the normalization condition J w(X, pL, u) dX = 
1. It is worthy to note that in view of homogenity of Dirac delta- function the tomogram has 
the same homogenity property. Connection between the density matrix p(x, x') and quantum 
tomogram w(X, fi, u) may be expressed through the following relation 
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4 Evolution equation in tomographic representation 



In this Section we introduce the quantum kinetic equations in tomographic form. The general 
evolution equation for the density matrix of the open system was suggested in P^ITdIITo] !^]. 
Solution of this equation corresponds to the Hermitian nonnegative operator with constant trace 
at all moments of time. According to ^H] the most general evolution equation which keeps stated 
properties takes the form: 
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where H is Hamiltonian, Vj are arbitrary linear operators and n is an arbitrary integer. This 
equation, presented in coordinate representation, has the form: 
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(35) 



where p = —id/dx, p' = —id/dx'. The kinetic equation in such a form can be easily transformed 
into equation on quantum tomogram. Indeed, the density matrix in coordinate representation can 
be expressed through tomogram w(X,p,, v) using relation ([3*3^1. Therefore we can establish the 
following correspondence of the action of the operators on the density matrix p (x, x') and the 
marginal distribution w(X, /x, v) [3P| : 
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One can see that all the operators in right hand side of above formulas have scale invariant, i.e. 
they are not changed if we replace Xj — > AjJQ, /ij — > Aj/ij, z/j — > AjZ/j. Thus using (|33j) and ([3*B|l one 
gets what we call evolution equation in tomographic (or probability) representation 



-i (H (Q, P)-H (Q', -P')) + 2 E ( 2 ^' W« P ) V i W« - p/ ) " 



V$ (Q, -P) (Q, P) - Vj (£?', -P') V 3 (Q', -P') 
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w(X,^,v,t), 



where operators Q, Q , P and P act on ^t, v) according to (JHEJ). Due to scale invariance of 

(|36|) the evolution equation preserves the homogenity property of tomogram. 
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5 Partial case of one-dimensional evolution equation in 
probability representation 

There exist several examples of kinetic equation for open systems. In this section we consider the 
tomographic evolution equation for the case of one- dimensional oscillator with Hamiltonian H of 
the form 

H = p l + - (38) 



(we assume that H = 1, m = 1, uj = 1) and the only operator V(q,p) = ux + vft, taken as a 
linear combination of the operators x and ft with arbitrary complex coefficients u and v. Some 
aspects of tomographic representation for this oscillator were considered in [18]. For u = \f^/2 
and v = i a/7/2 we have the standard quantum kinetic equation 
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For such a choice of operators H and V(q,ft) equation (}3T)) for the tomogram w(X,fi, u,t) can be 
represented in the following way: 
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w(X, /i, v, t). 



To solve this equation we make the Fourier transformation 
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and obtain the following equation for Fourier components w(k, p,, u, t) of the tomogram w(X, p,, u, t): 



w(k, p,, v, t) = — i ( ^ r£ — i Im (uv*) ) w(k, p,, v, t), 



(42) 



where 4- vector £ = (p^,p v ,x^,x u ) = (—id/dp,,—id/dv,p,,v) and 4x4 symmetric matrix T is of 
the form 



r 



r r 

J- pp L px 

r r 

1 xp 1 XX 



-lm(uv*) 1 

-1 -Im(ui;*) 
— lm(uv*) —1 —i\v\ 2 k 2 iRe(uv*)k 2 

1 — lm(uv*) iRe(uv*)k 2 — i\v\ 2 k 2 



(43) 



Equation (|4^j) for the Fourier components w(k, p,, u,t) is Shrodinger type equation with effective 
Hamiltonian given by a quadratic form. Equations of such a type can be solved in the framework of 
time dependent invariants method [31j. According to this method, Green function G(p,, f,pf, v' ,t) 
of equation i.e. function which connects w(k, p,, u, t) and w(k, p,, u, 0) by the equality 



w(k, /x, u, t) — / G(/i, z/, /i', z/', t)w(k 1 p! ', z/', 0) d// o?z^', 



(44) 
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can be obtained by solving the system of linear differential equations on 2 x 2 matrices Ai, A2, A3, 

Ai = AiT xp — A 2 r pp , A 3 = A 3 T xp — A^Tpp, 

(45) 

A2 = AiT^ — A 2 r px , A4 = A^T XX — A4r px , 

with initial conditions Ai = E 2 , A 2 = 0, A 3 = 0, A 4 = E 2 , where E 2 is 2 x 2 unity matrix. The 
solution of this system for the matrix T given by (|43|) reads 
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cost 
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a b 
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a — c 






sinh^— Im(-uf *)t \ cost+ 

cosh^— lm(uv*)t^ sint+ 
sinh ( — Im(-uf *)t ) sin t, 



(46) 



where constants a, b and c are determined by 
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For the case Im(wt> *) = the constant a and c are not defined but in the expression for matrix 
A 2 remains the only term in which these constants appear in the combinations a + c and a — c. 
It is easy to see that these combinations are defined even if lm(uv*) = 0. The Green function 
G(/i, v, //, z/, t) of equation (j4^|) for the case of A-matrices (f4T)|) is determined by formula pH] : 



1/, //, z/',t) 



S[ (fi, u) - A^/i', v') ) exp --(//, u')A 2 A 1 



(48) 



Vdet A 4 

Let us consider the evolution of the coherent state of the system with Hamiltonian ()38|) . Notice, 
that obtained Green function allows consideration of evolution of an arbitrary state. The tomogram 
corresponding to the coherent state of the oscillator is (sec e.g. 32J) 



w a (X, n, v, 0) 
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V 



yU 2 + Z/ 2 



(49) 



Taking the Fourier components of the tomogram w a (X, /i, u, t) one can find the evolution of 
w a (k,fi, v,t). After performing inverse Fourier transform of w Q (k,fi, u,t) we obtain the evolution 
of the tomogram w a (X, fi, u, t) of the form: 

[X - A// - Su] 2 



w a (X,fi, v, t) 
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y/ir(Cfi 2 + Dv 2 + Efiu] 
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: exp 



C/z 2 + Dp 2 + £/zz/ 



(50) 



where X = (x), 5 = (p), C = a xx , D = a vv and E = 2a xp . Considering evolution of the 
coherent state, one should distinguish two cases: lm(uv*) ^ and Im(w*) = 0. The coefficients 
C — C(u,v,t), D — D(u, v, t) and E = E(u, v, t) reads 

For the case of lm(uv *) 7^ For the case of Im(uv *) = 

C = (l -d + ccos2t - esin2t)x C = 1 + (\u\ 2 + \v\ 2 )t- 

x exp (2\m(uv*)t^ + d - c, - ^(l«| 2 - M 2 ) sin2t, 

D = (l -d-ccos2t + esin2tJx D = 1 + (|w| 2 + \v\ 2 )t+ ^ 

x exp (2lm(uv*)t) + d + c, + ^{\u\ 2 - \v\ 2 ) sin2t, 

E = -2^csin2t + ecos2t) x E = 2(\u\ 2 - \v\ 2 ) sin 2 t, 

x exp (2lm(uv*)tj + 2e, 

where constants c = c(u, v ), d = d(u, v) and e = e(u, v) are defined as 

(\u\ 2 — \v | 2 ) Im(w*) — 2 Re(uv*) 



2(1 + lm 2 (uv* 



\u\ 2 + \v\ 2 
2 lm(uv*) ' 



(52) 
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2^1 + Im 2 (W 

For both cases the coefficients A = X(u, v,t) and 5 = 5(u,v,t) can be represented as 

A = v^^Re (a) cost + Im (a) sinij exp (lm(uv*)t\ , 
5 = V2 ^Im (a) cost — Re (a) sintj exp ^Im(-uf *)t\ . 

Let us find the parameter /x , that characterizes the purity of the state w a (X, /1, z/), or the so called 
purity parameter. In terms of tomograms it can be written as [HO] : 

1 f 

Ho = — / w a (X, /J,, v, t)w a (Y, — //, — v, t) exp {i(X + Y)} dX dY d\i dv. (54) 
2n J 

Evaluating this integral one obtains the following expression for the purity parameter: 

1 . . 

it = ; ■ (55) 

^CD-E 2 /A 

Considering this expression one concludes that for parameters u = y/j and v = iypy the purity 
parameter equals to 1, i.e. for such a choice of u and v the ground and coherent states of the 
oscillator remain pure for all moments of time. For another choice of these parameters the purity 
/j,q is not equal to 1. For example the dependence of the purity on time for some other values of 
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(a) u = 1, v = lOi 



(b) u = 1, v = 1 + 2i 



Figure 1: Dependence of the purity on time for particular values of parameters u and v 



parameters u and v is shown on the Figure H Analyzing the behavior of fio depending on the 
values of fi and v we are able to make some general conclusions. If lm(uv*) = or lm(uv *) > 
then the purity /xo tends to the when the time t tends to the +oo. If lm(uv*) < then the purity 
Hq tends to the limit 



fx (u,v, +oo) 



a/c? 2 (m, v) — c 2 (u, v) — e 2 (w, v ) \ 



4 Im (uv*) 1 + 4 Im 2 (mw* 



(M 



If I 2 ) + 16 Im (-uf * 



(56) 



From the inequality Im 2 (-ut>*) ^ \u\ 2 \v \ 2 it is obvious that this limit is always less than unit and it 
is equal to the unit if | lm(uv*)\ = \u\\v\ and \u\ = \v\. 



6 Conclusion 

To summarize we point out the main results of this work. We reviewed the generic evolution equa- 
tion for density operator of open system in probability representation of quantum mechanics. In 
this representation the state evolution is described by conventional probability density (tomogram) 
evolution. In this sense the obtained equation is a generalization of Moyal evolution equation |22] , 
written for Wigner function • In our case we set the generic quantum kinetic equation for open 
systems similar to the Fokker-Plank type equation. The obtained equation is written in terms of 
quantum state tomograms. 

As example of such generic equations we have shown that the damped oscillator evolution 
depends on the character of chosen values of parameters available in the collision terms of the 
kinetic equation for the oscillator. 
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